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The unique capabilities of capacitance measurements in bilayer graphene enable probing of layer- 
specific properties that are normally out of reach in transport measurements. Furthermore, ca- 
pacitance measurements in the top-gate and penetration field geometries are sensitive to different 
physical quantities: the penetration field capacitance probes the two layers equally, whereas the 
top gate capacitance preferentially samples the near layer, resulting in the "near-layer capacitance 
enhancement" effect observed in recent top-gate capacitance measurements. We present a detailed 
theoretical description of this effect and show that capacitance can be used to determine the equilib- 
rium layer polarization, a potentially useful tool in the study of broken symmetry states in graphene. 



I. INTRODUCTION 



Capacitance measurements probe the energy cost of 
moving charge between different parts of a system. In 
a classical system, this energy cost is a purely geomet- 
ric quantity and consists of the electrostatic energy. In 
contrast, capacitance measurements performed on quan- 
tum systems can access a range of subtle and interesting 
phenomena. In particular, Pauli exclusion in degener- 
ate electronic systems gives rise to a characteristic quan- 
tum contribution to the internal energy. The associated 
contribution to capacitance, known as 'quantum capaci- 
tance—, is proportional to the electronic compressibility 
^j. In addition, at low carrier densities, the internal 
energy is dominated by electronic correlations, result- 
ing in a so-called negative compressibility contribution 
to capacitance^. In low dimensional systems these ef- 
fects can amount to a sizeable contribution, making ca- 
pacitance measurements a powerful probe of many-body 
effects^. Moreover, whereas electrical transport is often 
dominated by a small subset of electronic states, capac- 
itance probes all states equally. Consequently, capaci- 
tance is a useful tool in the study of phenomena in which 
localization plays a role, such as quantum Hall effects and 
the metal-insulator transition^—. Under certain condi- 
tions, the quantum capacitance can become an order-one 
effect. 

Graphene and its bilayer are ideal materials for the 
application of the capacitance technique. The two- 
dimensional geometry of these materials permits the 
placement of proximal metal gates^— , electrolytic 
solutions^, or scanning probe head s 14 i 15 , all of which 
can be used to probe capacitance. Interesting results 
have been obtained for monolayer graphene, in which 
this quantum capacitance was found to dominate the to- 
tal capacitance near the Dirac point even at room tem- 
perature^. Surprisingly, the compressibility measured 
at low temperature^ was found to be well described 
by the noninteracting massless Dirac model, a fact at- 
tributed to an exact cancelation of correlation effects in 
the monolayer—. In bilayer graphene (BLG), in contrast, 
the interaction effects are expected to be strong, po- 



tentially leading to novel many-body states near charge 
neutrality-^— . Such effects, if they exist, would directly 
manifest themselves in compressibility measurements 2 ^ 

In this paper we discuss the unique capabilities of ca- 
pacitance measurements in BLG. Due to the finite in- 
terlayer separation, capacitance measurements can probe 
layer-specific properties that are out of reach in conven- 
tional transport measurements in which the layers are not 
contacted separately. Motivated by recent experiments, 
we calculate the effect of a gate-induced charge imbal- 
ance between the layers on the measured capacitance in 
several geometries, taking into account electronic inter- 
actions and short range disorder. We interpret the pe- 
culiar electron-hole asymmetry observed in top-gate ca- 
pacitance measurements^ in terms of a "near-layer ca- 
pacitance enhancement" , which is a combined effect of 
van Hove singularities (vHs) in the BLG band structure 
and the interlayer screening. We show that capacitance 
experiments can be used as a which-layer probe, offering 
a unique capability in studying electronic properties of 
graphene. 



II. THE NEAR-LAYER CAPACITANCE 
ENHANCEMENT 

Recently, capacitance techniques have been applied to 
dual-gated bilayer graphen o 11 ' 12 . The geometry of these 
devices allows the electrostatic potentials on the two lay- 
ers to be varied independently, enabling independent con- 
trol of both carrier density and the gap in the electronic 
spectrum 2Zr2&. In the absence of external fields, BLG 
is a metal characterized (at sufficiently low energies) by 
approximately parabolic valence and conduction bands 
which touch at the corners of the hexagonal Brillouin 
zone (at the K and K' points). The degeneracy at this 
band crossing is protected by the symmetry of the BLG 
crystal structure, in which atomic sites on different layers 
are equivalent under transformations of the point sym- 
metry group. Application of an external electric field per- 
pendicular to the layers breaks the which-layer symmetry, 
turning BLG into a semiconductor with a gate-tunable 
band gap. At not too strong fields the gapped state can 
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be described^ by projecting the tight binding Hamilto- 
nian on the low-energy subspace of wavefunctions (i/'i, 
1^2) where the subscript indicates the layer index, giving 
the two-band Hamiltonian 
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P± = Px ±iPy, (1) 



where momentum p is measured relative to the K (or 
K') point and V\, V2 are the potentials on each layer, 
controlled by external gates or dopants. The Hamiltonian 
(fTj) features a band gap of size A = \v\ — U2I, and a pair 
of vHs in the density of states of inverse square root form 
positioned on cither side of the gap at e = vi and e = v 2 - 
The field-induced gapped state is characterized by in- 
tcrlayer density imbalance, in which the occupancies of 
the two layers are very different for v\ = v 2 and for 
v\ 7^ V2- For the balanced bilayer (v\ = 1*2) the wavefunc- 
tion amplitudes on each layer are equal (up to a phase); 
however, in the presence of an imbalance (v± 7^ v%) the 
amplitudes become unequal. This leads to population 
imbalance between the two layers, 
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with a higher occupancy on the layer which has lower 
energy. This layer population asymmetry results in a 
strong asymmetry in the partial (layer specific) densities 
of states: since each vHs shows up only in the partial 
density of states for one of the two layers, the correspond- 
ing divergent contribution to compressibility comes only 
from the vHs-bearing layer, remaining finite for the other 
layer. 

As we discuss in detail below, the layer population 
asymmetry, Eq.Q, manifests itself in capacitance mea- 
surements. This is illustrated in Figfljb), in which top- 
gate capacitance found using a self-consistent model (see 
Sec lIV|) is plotted as a function of gate voltages Vt and 
Vb- The enhancement in capacitance associated with the 
band edge is stronger when the divergent vHs-bearing 
layer is facing the gate used to measure capacitance (top 
layer for C t and bottom layer for Cb in Fig. Q]a). We refer 
to this behavior as 'near-layer capacitance enhancement' 
(NLCE). This NLCE effect is seen in the capacitance 
map shown in FigQJb): the dark region, corresponds to 
the insulating state realized when the chemical poten- 
tial is positioned inside field-induced gap, is bordered on 
one side by a bright fringe corresponding to the NLCE. 
The markedly different contrast between the van Hove 
singularity- associated features positioned on cither side 
of the dark region, is associated with the density piling 
up on the near layer rather than the far one. 

This behavior explains the asymmetry observed in top- 
gate capacitance measurements [l2j ]. in which a feature 
identified with the vHs was observed only for electrons 
(holes) when the high (low) energy layer was nearest 
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FIG. 1: a) Bilayer graphene capacitor schematic. Layer den- 
sities (ni and 712) and electrostatic potentials (vi and V2) are 
controlled by voltages on external gates (v% and fb), which 
couple to the bilayer through the fixed geometric capacitances 
Ct and C°. Capacitance measurements^ are performed by 
measuring the current flowing through both layers in the pres- 
ence of an AC driving potential on one of the gates, b) Top 
gate capacitance as a function of external gate potentials for 
a clean bilayer, calculated using the self-consistent approach 
of Sec|lV][see Eq.® as well as Eqs.([T7 )l -ipO ]l and iP ]) -|[T5 ]) ]. 
The capacitance, which is small in the insulating regime and 
high in the metallic regime, is enhanced at the edges of the 
metallic region due to the presence of van Hove singularities 
in the density of states at the band edge. The enhancement 
is asymmetric, reflecting the asymmetric population of the 
layers, Eq.©. 

the gate from which capacitance was measured. In con- 
trast, no such asymmetry is expected for the capacitance 
measured using 'penetration field' geometry^, because 
the penetration field capacitance is more symmetric than 
the one-sided (top or bottom) gate capacitance. Indeed, 
no NLCE-type asymmetry was observed in the measure- 
ments reported in Ref . |LU] . As we shall see, the gate ca- 
pacitance and the penetration field capacitance measure 
fundamentally different characteristics of the system. Si- 
multaneous measurements of gate and penetration field 
capacitances can thus provide detailed and direct infor- 
mation on layer polarization of the bilayer. 

The NLCE effect is sensitive to the form of the vHs, 
which depends on the specifics of the dispersion relation. 
The simplest model for BLG, which we focus on below, is 
that of quartic dispersion, described by the Hamiltonian 
(HJ). A more detailed analysis^— , based on the four 
band model, leads to a 'Mexican hat' structure in band 
dispersion near points K and K' . However, the Mexican 
hat dispersion and the quartic dispersion both lead to 
an inverse square-root vHs at the band edge, resulting in 
essentially identical NLCE effects. 

In this paper we develop theory of the NLCE effect. In 
section III wc calculate, using a two band model of BLG, 
layer-indexed densities of states, i/y — ~drii/dvj, where 
i,j,= 1,2 refer to the two layers. In section IV we de- 
velop a many-body approach that describes interactions 
of particles in BLG with other particles and also with 
gate potentials. Using a self-consistent Hartree-type ap- 
proximation, we derive expressions for several quantities 
of interest relevant to capacitance measurements in terms 
of the matrix elements Vij . Wc find that different cxper- 
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imental observables exhibit very different behavior. In 
particular, the gate capacitance exhibits strong particle- 
hole asymmetry and the NLCE effect (sec Fig.l), while 
the penetration-field capacitance is nearly particle-hole 
symmetric. In section V, we consider the effect of disor- 
der, and show that the asymmetry persists for relatively 
high disorder concentrations corresponding to the exper- 
imental regime. Finally, wc conclude with a discussion of 
the usefulness of different capacitance measurements in 
bilayer graphene for probing the layer-pseudospin texture 
of possible broken symmetry phases. 



III. THE VAN HOVE SINGULARITIES AND 
COMPRESSIBILITY IN CLEAN BLG 

The main features of the compressibility of BLG in an 
external field can be understood in terms of the many- 
body Hamiltonian 



(3) 



where H is the single-particle Hamiltonian ([1} and sum- 
mation over four flavors a = 1,2,3,4 accounts for the 
spin and valley (K, K') degrees of freedom. The inter- 
action is written in terms of density harmonics on the 
layers, n i)k = E P , a ^p,a^p+k (* = !> 2 )> 
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with Vk and Vk the intralayer and interlayer Coulomb 
interaction, 
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where d f=s 0.3 nm is the interlayer spacing in BLG. 

We analyze quantum corrections to the capacitance 
of gated BLG described by the Hamiltonian ^ us- 
ing a Hartree-type approximation. This is done in two 
steps. We first find the compressibility matrix of non- 
interacting fermions, formally setting ifi nt = in Eq.Q. 
In doing this, the BLG potentials v% and «2 are treated 
as external parameters. Next, in Sec lIVl we restore the 
interaction -ffi n t, adding to it the interaction between all 
charges, including those on the gates. We relate poten- 
tials i' 1 ( 2 ) to charges on the gates and the graphene bi- 
layer self-consistcntly, and use these relations to evaluate 
capacitance as a function of external gate voltages. 

The Hartree-type analysis presented in this paper does 
not account for correlation effects; however, estimates of 
the correlation energy and the analysis of compressibil- 
ity of BLG presented in Ref.[25[ indicate that the cor- 
responding correction to capacitance is small, except at 
very low values of disorder and temperature, where the 
BLG system develops an instability towards a correlated 
state. 



In recent experiment s 1 1 > 12 electronic states with dif- 
ferent doping relative to the neutrality point are probed 
by varying the potentials v\ and V2 through their re- 
sponse to the potentials Vt and applied to external 
gates. Metallic and insulating conductance regimes oc- 
cur when the Fermi level lies inside or outside the gate- 
induced gap22r— . The insulating regime was observed to 
accompany a drop in compressibility. 

It is convenient to introduce layer-symmetrized poten- 
tials v± = ^(fi±t>2). Within the two band model ([T]), the 
gap size is A = 2|w_| and the position of the gap center 
relative to the Fermi level is w+ — fi; the metallic and in- 
sulating regimes in a clean bilayer are then described by 
\v+— fi\ > \v-\ and /i| < \v-\, respectively. In exper- 
iment a 11 ^ 12 capacitance was measured with the graphene 
bilayer grounded. This situation can be described by a 
Fermi level pinned to zero energy, (i = 0. 

Particle densities on the two layers can be expressed 
as sums over all occupied states, 
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where /(p) = l/(e^ + 1). In what follows, we focus 
on the case of zero temperature, /(p) = 9(— e(p))- Using 
the eigenstates of the Hamiltonian ([I} and defining layer- 
symmetrized densities n± = ri\ ± ri2, we find 
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where A is an ultraviolet cutoff of order the band- 
width. Here vq = 2me 2 / (irti 2 ) accounts for the four-fold 
spin/ valley degeneracy, and can be written as 2 /was, 
where as is the Bohr's radius of BLG. The two cases 
in Eqs.©,©, metallic and insulating, correspond to the 
regimes > and < 

Using these expressions we can compute the entries 
of the compressibility matrix = —drii/dvj. The ex- 
pressions have different form for \v + \ > and for 
\v+\ < \v-\: 
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where we defined 
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with e = 2.71828.... Expressions ©-([III) are plotted in 
the left panel of Fig. [5] Note that the compressibility 
matrix is symmetric, = v |_. 

Different elements of matrix v have different physi- 
cal meanings. The diagonal element z/ ++ = —dn+/dv+ 
is the total charge compressibility. The diagonal ele- 
ment v = —dri-/dv- is layer polarizability. The off- 
diagonal elements v |_ = = —dn-/dv+ describe the 

charge-flavor response. The latter quantities are particu- 
larly useful, as they measure the layer distribution of in- 
cremental additions of charge, giving information about 
the layer polarization of the ground state: the quantities 

v |_ and are zero for an unpolarized bilayer, but 

nonzero in the presence of a charge imbalance. 

Rewriting Eqs.([7]),([5|) in terms of variables character- 
izing individual layers, rii, ri2, we obtain 
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Expressions (|13[) - (|15[) arc invariant under simultaneous 
1h2 exchange and gap inversion, r>_ — > — V—. 

Both of the diagonal compressibility matrix elements 
(fn and V22) exhibit an inverse square root divergence 
at the charge gap edge, where the density of single par- 
ticle states has a van Hove singularity. The two diago- 
nal compressibilities behave asymmetrically, diverging on 
opposite sides of the gap: dn\/dv\ diverges at v\ —t 0, 
while 8122/ dv2 diverges at v 2 0. In contrast, the off- 
diagonal compressibilities (i 7^ j) remain finite on either 
side of the charge gap and are symmetric (see Fig. [5J left 
panel). Inside the charge gap, |u + | < the diagonal 
and off-diagonal compressibilities are constant: 
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exhibiting no divergence at the gap edge. 
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FIG. 2: Energy dependence of the inter layer compressibility 
matrix elements Vij in the 1/2 (left panel, Eqs. (|13|) - (|15[ )) 
and v± (right panel, Eqs. (f9|l- (|ll|l ') bases for fixed interlayer 
asymmetry v~ = 50meV and A = 5eV. In the left panel, 
single layer charge compressibilities v\\ and V21 are divergent 
only on one side of the charge gap, allowing the interlayer 
asymmetry to be probed by single side capacitance measure- 
ments. In the +/— basis, this asymmetry is reflected by the 
charge-flavor response, 



gates, characterized by potentials v t and «b, charge densi- 
ties tit and rib, and geometric capacitances to the bilayer 
C t ° and C°. The bilayer is described by the potentials 
Vi and V2 and charge densities m and 72,2 induced by 
the external gates on the individual layers. Electrostatic 
energy of the bilayer itself is taken into account by in- 
cluding an interlayer capacitance CblGj which can be 
estimated from the "geometric" value obtained for a par- 
allel plate capacitor, Cblg = (47rd) _1 , with d ~ 0.3 nm. 
This electrostatic model amounts to the approximation 
that the charge density on the bilayer is of the for m 
n(z) — mS(z — d/2) + n 2 S(z + d/2). While corrections 
are expected due to the finite extent of the wavefunc- 
tions, these corrections amount, for the most part, to a 
renormalization of CblGj upon which our results do not 
sensitively depend. 

The quantities of interest obey the general electrostatic 
charge field relations 
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C°(v2-v h ) = -(n t +ni+n 2 -n h ), (19) 



n t + rii + n 2 + n h = 0. 



(20) 



IV. SELF-CONSISTENT CAPACITANCE 
CALCULATION 

We shall focus on the geometry pictured in FigQJi, 
which describes a dual-gated graphene device of the type 
studied in Rcfs.[l2[ and [ll|. The experimental system 
consists of a bilayer graphene sheet placed between two 



To complete the system of equations for charge densities 
and potentials, a set of constitutive relations for BLG 
must be used. These relations, which are of general form 
m = U2)) ri2 = f 2 {vi 1 V2) 1 will be calculated in 

subsequent sections. 

Capacitance measurements are done in the finite fre- 
quency regime, by applying a small AC bias (on top of 
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the DC bias used to control density and intcrlayer imbal- 
ance) to one terminal of the device and then recording the 
resulting change in charge density on a second terminal. 
Choice of terminals distinguishes top (back) gate capac- 
from penetration field capacitance, C p , 
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After eliminating n t and rib from Eqs. (fT7)) -(f2"0 |) by ex- 
pressing them in terms of other variables, n t = C®(vt — 
vi), rib = Cb( w b — V2), the remaining two equations are 
linearized with the help of the matrix of inter- and in- 
tralayer compressibilities 
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FIG. 3: Calculated Ct (a) and C p (b) for the clean bi- 
layer. Different color traces correspond to different values 
of the top gate capacitance, measured relative to a fixed 
C° (taken to be 120 aF/^tm 2 corresponding to the stan- 
dard 285 nm S1O2). Penetration field traces are normal- 
ized by the geometric value corresponding to full penetration, 

c% = {i/c h , + i/c^ G + i/c?y 1 . 



where C is a matrix of geometric capacitances, 
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These expressions account for both the geometric and 
'intrinsic' capacitance of BLG. 

Solving for 5vi, Sv2, we find the charges induced on 
each layer by the gate potentials: 
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Here the first term describes the geometric capacitance, 
which would be the only contribution if the electronic 
system in BLG was infinitely compressible, v — > 00. The 
term proportional to —C{v + C) -1 describes the quan- 
tum capacitance contribution. Combining equation (|25[) 
with the relations for rit and rib, all three capacitance 
observables can be calculated: 
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(28) 



These quantities implicitly depend on the gate potentials 
through the compressibility matrix Vij . 

Notably, different capacitance observables depend on 
different combinations of the compressibility matrix ele- 
ments, and obey different symmetries. The penetration 
field capacitance C p is dominated by the off diagonal 



component of the (necessarily symmetric) compressibility 
matrix. As a result, for a symmetric device (C^ = C°) it 
is invariant under interchanging layers 1 and 2 and there- 
fore does not exhibit the NLCE effect. In contrast, the 
expressions for Cb and Ct are not 1 <-> 2 invariant. In 
particular, the last term in the expression for Ct, propor- 
tional to V22, changes to v\\ upon layer permutation. As 
shown in the previous section, in the presence of a layer 
imbalance these two quantities are not the same, leading 
to the observed NLCE observed in Reffl2l. 

In a device in which all capacitances can be measured, 
combinations of the measured quantities can be combined 
to probe the charge-flavor response. For the simplest case 
of a symmetric gate configuration (£7° = C t °), 
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+ H 



Because this quantity is proportional to v |_, it can 

be used to probe both gate-induced and sponta- 
neous layer polarization, allowing direct experimen- 
tal measurement — somewhat analogous to Knight Shift 
measurements for spin — of the ground state layer polar- 
ization. 



V. THE EFFECT OF DISORDER 

In the devices used for capacitance measurements in 
Refs.[12|.|ll|. graphene flakes were supported by a silica 
substrate. The carrier mobility in such devices was of 
order 1,000 cm 2 /V sec. For such low-mobility devices, 
taking into account the effect of disorder is crucial for 
developing a sensible model of the experimental data. 
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FIG. 4: The effect of disorder on the density of states. Partial 
density of states p; , Eq. (|38[) for layers i = 1 (solid lines) and 
i = 2 (dashed lines) of a graphene bilayer, obtained from the 
self-consistent Born approximation, Eqs. (l36|) , ll37[) . Increasing 
the disorder strength leads to smearing of van Hove singular- 
ities and, eventually a closing of the energy gap. 

Full quantitative description of experiments requires in- 
cluding realistic disorder, which is likely long ranged—, 
along with the effects of electronic correlations^ which 
can give quantitative corrections to the electronic com- 
pressibility. However, the the key features of the data 
are captured by a simpler short range disorder model^l, 
which involves delta-function impurities localized on car- 
bon sites: 



(30) 



with potential u(x) = J^i U8(yi — Xj) taking values U on 
the carbon sites occupied by impurities, and zero else- 
where. The impurities are assumed to be distributed 
randomly with concentration n. 

The problem ([50]) can be analyzed using a self- 
consistent T-matrix approximation (SCTA). The SCTA 
approach provides a somewhat more general approach 
than the self-consistent Born approximation, and is re- 
duced to the latter for weak disorder. 

We evaluate the DOS and the total energy by em- 
ploying disorder-averaged Greens functions expressed 
through the layer-indexed disorder-averaged self-energies 



G(e,p) 



Vi - Si 



-*p 

v 2 - S 2 



(31) 



where t p is the kinetic energy operator- 
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e -™=i _|_ e~"*"" 2 ) 2 . An infinitesimal imaginary part ±i0 
should be added to e to obtain the retarded and advanced 
Greens functions. 



The self-energy is approximated by the average values 
of the T-matrix, evaluated separately for the sites on 
layers 1 and 2, 



Ei(e) = ft(T 1 (e)) ) E 2 (e) = n(T 2 (e)). 



(32) 



Here n = npo is the adatom density with po = 2/3 V3a 2 
the density of type 1 sites. The quantities Ti( 2 ), written 
as a 2 x 2 matrix, are given by 
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where U = U/pq. For realistic values of v\ and v 2 the 
integral of the Greens function over the Brillouin zone is 
dominated by the regions near K and K'; approximating 
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(34) 

This expression is valid for £1(2) small compared to the 
bandwidth. Combining this result with Eq. (|32[) , we ob- 
tain two coupled equations for E\, e 2 : 



£l = £ — V\- 
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(35) 



where we defined A(e) = ^/^i/ £ 2 and P = mU/2. Solving 
these equations for E\, e 2 as a function of e, we find the 
Greens function (|3"Tj) and use it to calculate the density 
of states, 
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where the integral is identical to the one in Eg. (|3~4"1) . A 
factor of two was inserted after integration to account for 
spin degeneracy. 

The density of states is expressed through the quantity 
A(e). Taking the ratio of the self-consistent equations 
for Ei and e 2 , Eq. p5p . we obtain a single equation for 
the quantity A. Focusing on the case of weak disorder 
potential and expanding in U, we arrive at 



A 



e — Vi + i"f/X 
£ — v 2 + 



7 



mil 
~2pJ 



(37) 



where the terms linear in U have been incorporated in 
the quantities vu 2 y Once X(e) is found from Eq. (f3"T|) , it 
can be plugged into Eq. (|3"6f to obtain partial densities of 
states on each of the layers (see FigJ3|), 



pi{s) = yRcA 1 , 



p 2 ( £ ) = ^ReA. (38) 



In the absence of disorder, 7 = 0, we have A = 
y/ (e — vi)/(e — v 2 ), which gives van Hove singularities of 
an inverse square root form at the band edges e = vi,v 2 
as found in section I. In the presence of disorder, these 
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c/c' c/c° 

t t p p 




v + (meV) v (meV) 

FIG. 5: Top gate (left panel) and penetration field (right 
panel) capacitance for different values of the short-range dis- 
order parameter 7, here measured in meV. Interlayer asym- 
metry parameter v~ — 50meV and the cutoff A = 5eV. Geo- 
metric parameters are chosen to match experiment reported 
in Ref.52, C t °/C* = 30, C£ = 120aF//im 2 . Color scheme 
corresponds to varying values of 7 as in Fig. 0] 



singularities are washed out to varying degrees. As shown 
in Fig|31 this washing out proceeds by both reducing the 
height of the vHs peak and closing the gap. Crucially, the 
'off '-layer density of states at the energy of the 'on' layer 
vHs peak increases with disorder. This has the effect of 
increasing the screening effect of the 'off' layer when it 
lies closer to the gate used to measure capacitance, en- 
hancing the NLCE effect for disordered samples. 

To calculate experimental capacitances, Eqs. (|27|) - 
(|2"5)l . the partial densities of states are integrated numer- 
ically with respect to energy and then rcdiffcrcntiatcd 
with respect to the appropriate energy variable, v\ or 
v 2 - In Figure El the results for both top gate and pen- 
etration field capacitance for a device with electrostatic 
parameters resembling those in Reffl2l are plotted. The 
asymmetry of top gate capacitance survives disorder av- 
eraging, and indeed is enhanced. For intermediate values 
of disorder, electrons and holes display qualitatively dif- 
ferent behavior: the non-monotonic vHs feature survives 
for holes but is completely obliterated for electrons, as 
observed in Reffl2l 



sitivity to the interlayer imbalance arises despite the fact 
that the layers are not contacted separately: the relative 
proximity of the layers to the top- and bottom- gates, 
combined with the interlayer screening, allows capaci- 
tance measurements to access layer specific quantities. 
Gate capacitance measurements preferentially probe the 
nearer layer, leading to the NLCE effect as the near layer 
screens the far layer. Consequently, in the presence of a 
layer imbalance, top- and bottom- gate capacitance mea- 
surements will be different. This difference is the sig- 
nature of layer polarization, allowing its unambiguous 
experimental determination. 

Our analysis provides an explanation of recent top gate 
capacitance experiments on dual gate bilayer graphene 
structure s] 11 ' 12 . Since the degeneracy of the band crossing 
in the BLG spectrum at the K and K' points is linked to 
inversion symmetry, the gate-induced density imbalance 
and the opening of a band gap go hand in hand22~— . 
As we have shown, this imbalance can be probed di- 
rectly through NLCE measurements; to our knowledge, 
the NLCE- type asymmetry observed in Reffl2l is the first 
direct experimental evidence of layer imbalance in BLG. 

The possibility of probing layer polarization directly 
through capacitance measurements has implications be- 
yond the study of gate-induced gap opening. Recently, 
experimental sample quality has improved to the point of 
allowing the observation of a multitude of novel features 
likely associated with electronic correlation a 26 ' 38 " — . A 
large number of possible broken symmetry states, arising 
in the presence and in the absence of magnetic field, have 
been explored in the theoretical literature^— , including 
several mutually exclusive scenarios for the ordering at 
low densities and small electric and magnetic fields. The 
main open questions pertaining to these states have to do 
with identifying broken symmetries and determining the 
exact structure of the order parameter and excitations. 
Future NLCE measurements, by offering a direct method 
for determination of the layer polarization, will help to 
narrow down the possibilities for these new states. 
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